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Abstract : The aim of this paper is to find higher order geometrical corrections to the Einstein- 
Hilbert action that can lead to only second order equations of motion. The metric formalism is 
used, and static spherically symmetric and Friedmann-Lemaitre space-times are considered, in four 
dimensions. The FKWC-basis are introduced in order to consider all the possible invariant scalars, 
and both polynomial and non-polynomial gravities are investigated. 


1 Introduction 

Most of the equations of motion describing physical effects are second order, that is, we need to 
specify either an initial and a final position in space-time to describe the dynamics between them, 
or we need an initial position and velocity to describe how the system will evolve. Concerning 
General Relativity (GR), for which the gravitational field g^i/{x) is encoded into the geometry of 
space-time : 


ds^ = g^,j{x)dx^dx '', ( 1 ) 

The Einstein field equations, describing the dynamics of the geometry, are also second order ones. 

However, it is well known that two of the simpliest solutions of GR, the Schwarzschild metric 
and the Friedmann-Lemaitre one, suffer from the existence of singularities. When one is dealing 
with ordinary matter, this is a general fact. Furthermore this theory alone is not able to describe 
dark energy, even though the inclusion of a suitable cosmological constant is sufficient. But then, 
other problems arise, like the cosmological constant one [T] and the coincidence problem. Therefore, 
one can think about modifying the Einstein equations, in the hope to describe dark energy and to 
cure singularities. In order to do so, one can add higher order invariant scalars, like , to 

the Einstein-Hilbert action to have high energy corrections that could describe what really happen 
around the singularities [a 13]. With regard to dark energy issue, see [1013 Hi- 

Within an higher order modified gravity model, the equations of motion will no longer be second 
order ones : there will be more than two initial conditions to specify in order to find the dynamics, 
so to keep the physical sense of what is an equation of motion, it is needed to introduce new fields 
for which these additional initial conditions would apply, such that at the end, the theory would 
involve two dynamical fields, with second order equations of motion for both of them. By doing so, 
we face an important problem, that is the presence of Ostrogradsky instabilities (see, for example 
0) : the new field defined in this way can carry negative kinetic energy such that the Hamiltonian 
of the theory is not bounded from below and can reach arbitrarily negative energies, what would 
make this theory impossible to quantize in a satisfying way [4]- And there are no general rules to 
avoid this problem, although a well known class of modified gravity, equivalent to GR plus a scalar 
field, the f{R) one, might not suffer from this problem 0. 
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Moreover, with a new field involved in the dynamics of gravity, this last would not be a fully 
geometrical theory anymore, which is yet one of the most important implications of General Rela¬ 
tivity. Nevertheless, it is possible to find second order equations of motion from the addition into 
the Einstein-Hilbert action of higher order scalars m- In this way, the Ostrogradsky instability 
may be avoided and there are no additional field involved in the dynamics, so these corrections 
can be said to be ’’geometrical” ones. This kind of modifications are the Lovelock scalars, but it 
turns out that in four dimensions, the only higher order scalar made of contractions of curvature 
tensors (only) that leads to second order equations [TU] is the so-called Gauss-Bonnet invariant : 

= + ( 2 ) 

which is however a total derivative in four dimensions m, and then it does not contribute to the 
equation of motion: 

= ^. (3) 

This result is background independent, which means that if we want to find a second order 
correction for all possible metrics in four dimensions, then this unique term does not contribute to 
the dynamics. That is why, in order to find anyway significant corrections to General Relativity 
that could cure some of its problems, we will search for additional terms that will give second order 
equations for only some specific metrics : the most studied ones, that suffer from singularities, the 
FLRW space-time describing the large scale dynamics of the universe, and the static spherically 
symmetric space-time describing neutral non-rotating stars and black holes. 

We note however that our way to find second order corrections is not at all the only possible 
one. There are other formulations of GR than the metric one, where the equations of motion are 
found by varying the action with respect to the metric field only. In the spirit of gauge theories, 
one can also vary the action with respect to the connections and independently with respect to the 
metric. Then, it is possible to find second order corrections with no a priori background structures 

[13- 

In some sense, our approach is similar to Horndeski’s theory which is the most general one 
leading to second order equations of motion for gravity described by a metric coupled with 
a scalar field (j) and its first two derivatives |13] . This theory involves non linear higher order 
derivatives of the scalar field, like (□</>) ^, and yet leads to second order. Moreover, if all the matter 
fields are minimally coupled with the same metric </>), one can expect the equivalence 

principle to hold [14], which is also a fundamental feature of GR that one wants to keep. 

Briefly, the outline of the paper is the following. First we consider all the independent scalar 
invariants built from the metric field and its derivatives, for example of the form (DR) , and 
see if some linear combinations of them, or, in the spirit of |15j and m, if some roots of these 
combinations, could lead to second order differential equations for FLRW space-time and static 
spherically symmetric one. The basis of independent scalars that are needed have been presented 
in but for specific backgrounds, we will show that this basis may be reduced. Furthermore, 
we will start to exibit, order by order for FLRW, the existence of polynomial and non-polynomial 
gravity models that give second order equations and polynomial corrections to the Friedmann 
equation. Finally, we will investigate the static spherically symmetric space-times. 


2 Order 6 FKWC-basis 

The basis of all independent invariant geometrical scalars involving 2n derivatives of the metric 
are separated into different classes, depending on how many covariant derivatives act on curvature 
tensors. For order 6 (n = 3), the first class, that does not involve explicitly covariant derivatives, 
from Cl to Cs, is denoted by Rg 3 : these scalars are built with six derivatives of the metric and by 
the contraction of 3 curvature tensors. The two other classes, 7?.°2 qj and contain scalars 

that involve respectively, a curvature tensor contracted with two covariant derivatives acting on 
another curvature tensor (from Jjfi to and two covariant derivatives, each acting on one cur¬ 
vature tensor (from Jjfg to CCs) : 
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r £1 = R^’'^^Rc.l3apR^P,,, 

’ £2 - A apR vpR pa 

1 £3 = 

, £4 = RR^^'^^^Rp.cp = RT 

] C5 = Rf^-'^^Rpc.R.p 

n _ iDiiv D do: 

, Lq — nr ^ 

[ £7 = RR^^''Rp^ = RS , 

Cs = R^ 

( Jfi = RUR 

^2 = Rp.aR^'' 


^4 = 


, = ^aRp.V'Rf^^ 

[ ^7 = 

, ^8 = VaRV^R 

Recall that our aim is to see that if we consider all these scalars, there are quite natural modified 
gravity Lagrangian densities that we can expect to lead to second order equations of motion and 


that actually do. There are linear combinations of all the scalars of the basis, but also for example 
square-root of the 72.°^ ^j-class, or cubic-roots of the 72g 3 one, even if we are not going to study this 
last because we search here for high energy geometrical correction to the Einstein-Hilbert action. 
We write down this fact as : 

^ -f RVVR + WRVR )++{/!:{«»). 

Because inside a same class, the scalars have approximatively the same terms in their expansions, 
we can indeed expect to cancel higher order derivatives for some specihc combinations of them, 
and then to have second order equations of motion. 

Now, let us write down some definitions that allow to find relations between these scalars, 
coming from the fact that we are going to restrict our study to specific backgrounds, the Friedmann- 
Lemaitre-Robertson-Walker (FLRW) metric and the static spherically symmetric one, both in four 
dimensions. 

For both of them, there are between the scalars relations coming from the Lovelock theorem 
(that are not taken into account in Ref [IZ]), and also for FLRW, relations coming from the fact 
that this is a conformally invariant flat metric. All these relations are written in the first Appendix, 
and we need to express them to define the Weyl tensor, as : 

— R/ii/CK/3 RfifUdua T RufUQRuctQ “f 0 9fil39i>ct)R- (4) 

And the following rank 2 tensor that is null in four dimension because of the Lovelock theorem : 

Lij,u = ——gfiu£4 + Ra'i + ‘2-RR^v = 0, (5) 

where = R^na R^'^p and Pfiy = Ru-yR'^^. Indeed, if we vary the Lagrangian associated with 
the Gauss-Bonnet invariant with respect to the metric field we find : 

=^9 ( 6 ) 

But as we saw in equation (3), this Lagrangian can be written as a total derivative in four dimen¬ 
sions, which means that its contribution to the equations of motion is identically zero. 


3 Friedmann-Lemaitre Space-time 

3.1 Order 6 

We start with the flat Friedmann-Lemaitre cosmological metric describing the dynamics of the 
universe at very large scale, in the simpliest manner : 

ds'^ =—dt^ + a{t)‘^(dr'^ + (7) 

where a{t) is the scale factor and dQ^ = d9^ -|- sin^ 9 d(j)^ is the metric of the 2-sphere. This 
metric is conformally invariant to Minkowski space-time, and from the relations written in the first 
Appendix, we can choose the following reduced basis of all independent order 6 scalar invariants: 
(£4, Hgj-Cy, jSfs, jSfs). We note that there is one scalar less than for a general conformally 

invariant space-time coming from the particular metric of FLRW for which there is the additional 
relation: 

£1 = i(-£7 + 2£4). (8) 
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3.1.1 Linear Combination. correction. 


With this metric, we can right the most general order 6 linear combination of all the independent 
scalars : 


J = 


(viCi + ^ ((Jiivi,Xi) a(t)® + a2{vi,Xi) a{t) d{t) + a3{vi,Xi) a{t)‘^ a{tf' a{tY 


a{t) 

\ 

+ a4{vi,Xi) a{t)^ + cr5{vi,Xi) a{t)'^ a{t)^ + ae(vi,Xi) a{t)^ d{t) d{t) 

+ a7{vi,Xi) a{t)'^ + asivi,Xi) a(tf a{t)'^ + crg{vi,Xi) a{t)* d{t) , 


where the expressions of the aj in terms of (vi,Xi) are presented in the second Appendix. Setting 
all of them to zero allows to check that our list of scalars is a basis. We can then impose vi = vg = 
^3 = W 5 = Vs = 2^2 = 2^4 = 2^6 = = 0 to take into account the algebraic relations we have found. 

Moreover, in this section, we are interested in linear combinations of order 6 scalars that lead 
to second order equations of motion. Therefore, we can also consider equivalence relations (up to 
boundary terms) between the scalars, and there are three of them that remain after considering 
the previous algebraic relations : 

J d ‘^ x^/^^3 ^ J *^8 ’ J = - J d '^ x ^/^ £^3 


and 



d'^x^f^ { 2^8 + SA - l^Ce + £ 7 ) • 


We check that there is no other one by deriving the equations of motion for the scale factor 
considering the lagrangian L = a{t)^ J , that we substitute into the generalized Euler-Lagrange 
equation for third order lagrangians: 



( dL \ 

1+ 1 



^dL\ 

dt^' 

b5a(3)v 


^dd) 

dt ' 

y dd) 


Finally we only need to consider the combination : 


(9) 


J = V4C4 + VqCq + V7C7 + Xs^^Sj 


and after deriving the equation of motion and imposing a simultaneous cancellation of the higher 
order terms, we find that there is only one linear combination leading to second order equations : 

Ji = - 7 C 4 + 2{6£6 + £ 7 ) = 72H{t)^ (sH{t) + 2H{tf^ , (10) 

where H{t) = d(t)/a(t) is the Hubble parameter. We note that this linear combination only 
involves contractions of curvature tensors, which is the kind of corrections expected to follow from 
quantum held theory. Therefore, considering the following action that could represent a natural 
high energy geometrical correction to the Einstein-Hilbert action for FLRW space-time, 


( 1 


Si = j d^xy^ i? + - 7 + 12 + 2 ) 

where is the Lagrangian density for matter, one hnd the acceleration equation : 

3iJ (t)2 -h 2H{t) -mvH{tf {H{if + 2H{t)) = -Sirp, 


,( 11 ) 


( 12 ) 


where p is the cosmic pressure. Then the equation of the conservation of the energy, with p the 
cosmic energy density, 

■^ + 3i7(t)(p-bp) = 0, (13) 

gives the following modihed Friedmann equation : 

ZH{tf - 36 = Sirp. (14) 


One can solve it, and choose only the solutions that reduce to the standard equation when p is 
small. However, the only solution we are going to see in this work for FLRW spacetime is coming 
from a non-polynomial correction involving order 8 scalars, but the fact that Si is unique and 
second order could be a sufficient reason to study its cosmological solutions. 
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3.1.2 Non-polynomial gravity. correction. 

Now we want to consider all the order 6 linear combinations that are perfect squares, in order to 
consider non-polynomial corrections that are the square-roots of these squares. Rewrite the most 
general linear combination in terms of the Hubble parameter H (t) : 


J= ^ ViC^+ ^ =3iai{vi,Xi) H{t)^+ a 2 {vi,Xi) H{t)‘^H{t)+a 3 {vi,Xi) 

1=4,6,7 i=l,3,5,8 \ 

+ a4{vi,Xi) H{t)^ + '^5ivi,Xi) +aQ{vi,Xi) H{t)H{t)H{t) 

-\-a7{vi,Xi)H{if +aii{vi,Xi) +ag{vi,Xi) . 


To find squares, we use the following general procedure that is usefull for FLRW space-time, and 
necessary for spherical symmetry : Take the higher order perfect square In the expansion 

of our square, each term will be multiply by H{t), so the only terms that can enter inside are those 
Kiit) for which Ki{t)H{t) and Ki{t)Kj{t) exist in the expansion of order 6 scalars. Because of 
this, we need to impose the conditions : ctq = 0 , cts = 0 and CT 4 = 0 , what give xi = X 3 = 0 
and V 7 = —Vi — 5w 6/12. Therefore, CT 5 = 0 and there are only two possible forms of squares made 
of order 6 scalars : 

{vi£i +Xj£fj) = {5H{t)H{t) , 


And : {viCi + Xj£fj) = + 6 H{t)H{t)^ . 

i,j 

It means that all their square-roots can be decomposed in the basis 
Moreover, the general Lagrangian density : 


IT.' 


= ^Hit) + 5H(t)H(t) + jH{tf 


(15) 


leads to second order differential equations for all (^, d, 7), so we do not need to impose additional 
conditions on these coefficients and there are then only 3 independent second order corrections 
that we can hnd in this way. 

Now let us see what are the actual perfect squares that one can find. Solving the natural 
conditions for respectively the first and the second kind of perfect squares, = 4 ct 3 CTi , dg = 
ar = 0 and = 4 CT 3 dy , ai = CT 2 = 0, one find the following ones : 


J 2 = 2£4 + 12£6 - 7£7 = -72 4iJ^ -h 3HH 


J 3 = (^Ci - 12 Cq - C^a{ha + 18/3) -h 6 (a -K 3/3) 


(16) 


= -72 ( 3(a -h 4:(3)H{t)H{t) + {a + i(3)H{t) 


The last one being a general formula that gives perfect squares for any value of (a, /3). As we just 
saw, there are only 3 independent square-roots of these squares, so from the general expression J 3 , 
we can choose the two last to be the squares given by (a = 0, /3 = v^/O) and (a = 1, /3 = O): 


Jsy — — —36 




(17) 


And : J 3,2 = ( 6 . 5^5 + 5 ( 6£4 - - £ 7 )) = -72 (3H{t)H{t) + ‘ 


(18) 
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Indeed, one can check that the following combination, with Ci the signs inside de squares. 


^0 + vi eiy'-Js^i + V 2 e2-\/-J3,2 =6^3v^ (a + 4/3) i/q cq + 4i/i ei + 3 -^21/2 e.-^H{t)H{t) 

+6^-\/2 (a + 3/3) t'o cq + ei + V2v2 e 2 ^H{t), 

Vanishes for V 2 = —Q;z/oeo/e 2 and vi = — 3-\/2/3z/oeo/ei: and so the relation becomes explicitly, for 
all value of (a, /3) : 


{qU - ISA - /:7)a(5a + 18/3) + 6 (a + 3/3) (aJSfa + 

- ^^-(6^5 + 5(6£4 - 12/:6 - ^ 7 )) = 0 - 

Therefore, this general formula for perfect squares depends only on J 3 1 and y—as we 
said. 

We note here that an interesting property coming from the existence of an infinite number J^{a, /3) 
of perfect squares for which the square-root can be decomposed in a small basis is that it gives 
some non-linear algebraic relations between the scalars of the FKWC-Basis, that reduce it in a 
non-trivial way, and allow to have a very small number of independent corrections. Indeed, solving 
the previous equation for £4, we find : 

£4 = ^( - gjf’s + S.if’s + 108£6 + 9£7 - Y/£f 8 (l 8 =Sf’ 5 - 5 £f 8 )eoei) • (19) 

We can now calculate the equations of motion for the 3 Lagrangian densities '/—J 2 , y^—Jsp and 
\/—J 3 , 2 - First, one can check that the last one is in fact a topological term that does not bring any 
contribution to the equation of motion. Moreover, the first two lagrangian densities give the same 
equation of motion 5Aa(t)d{t) = 0 for the first one, and 18-\/2d(t)d(t) = 0 for the second one. It 
means that they are equal up to an invariant scalar T for which y/—gT is a total derivative, 

y-(2£4 + 12£6-7£7) = + T , ( 20 ) 

such that we can in fact consider a unique scalar (let us choose \J made of order 6 scalars that 
leads to non vanishing second order differential equations. We can note that T cannot be equal 
to y^—J 3_2 because \/—J 2 contains an terms. Therefore T could be found considering higher 
order derivatives scalars, and other perfect powers. In this case, powers of four for order 12 scalars 
for example. There is then another non-linear relation between order 6 scalars and (possibly) order 
12 ones. 

Finally, to recapitulate this part, we have found that it is natural to consider only one perfect 
square made of order 6 scalars : = V^RV^R. As a result, the action : 



52 




R -I- v^Z—V^RS/aR 



( 21 ) 


leads to a unique second order if^-correction to the Friedmann equation, as it is easy to see using 
the same reasoning as in the previous section. 


3.2 Order 8 

Now let us study the linear combination and squares made of order 8 scalars that lead to second 
order equations of motion. We do not copy all the FKWC basis for general metric, but we name 
the scalars according to their position in Ref m where this basis is fully written. The reduced 
FKWC basis for order 8 scalars in FLRW space-time is the following : 


6 



























'/Ci = i?4 , /Cio = , ICll=RR^^‘'‘^^R^\PR,,pp = R{\Cl-C 2 ) 

/Ci2=T2 , Mi = RD^R , TWa = , TWg = 

< TWg = V^ni?Vpi? , TWe = VpV^Vc,i?V'^i?'"“ , TWio = (Di?)^ , TWn = 

7Wi2 = , 7Wi 4 = VpV^i?„/3V^V'^i?“^ , TWig = , Mw = R^4 

, A^2o = SDR , A^33 = -R^s 


We also introduce the definitions JCg = R^^R°'^R^^^pR papa , -Adis = DR^^ORf^’' and 
A^i 6 = VpVi/Ro-^V^V^i?'^'' that will be usefull later for static spherically symmetric space-times. 


3.2.1 Linear Combination. iJ® correction. 

Consider the sum of all independent order 8 scalars for FLRW space-time : 


J = Y. + Z) 

Here, we follow exactly what we did for order 6 scalars. We derive the equation of motion associated 
with the previous sum, and see what conditions on (vi,Xj) cancel the equation, such that we 
find the 10 equivalence relations that exist between the scalars of the reduced basis. Therefore, 
we can consider only the following independent scalars with respect to the equation of motion, 
(/Cl,/Cio,/Cii,/Ci 2 , All, Atii, Ati 2 ) and the reduced sum : 

J= Y Vi]Ci+ Y ^3^3- 

i=l,10,ll,12 1 = 1,11,12 

We derive its associated equation of motion and see what a simultaneous cancellation of all the 
higher order terms implies for the coefficients {vi,Xj) : we find that the unique linear combination 
of order 8 scalars for FLRW space-time that leads to second order equation is : 

J 4 = /Cl - 48 /Cii - 9 /C 12 = 1728 ( iJ® -f 2 HH^ ). ( 22 ) 

Therefore, one may consider the action : 


^3 




R+13(^R^ 





(23) 


and see that it brings an i7® correction to the Friedmann equation. We note that this correction 
involves only contraction of curvature tensors like for the order 6 case. 


3.2.2 Non-polynomial gravity. correction. 

Correction to the Einstein-Hilbert action 

To hnd non-polynomial second order models from order 8 scalars, we follow exactly what we did 
for the order 6 and find the same kind of result : there are only two classes of perfect squares in this 
case. Those for which the square-roots give topological scalars, that does not give any contribution 
to the equation of motion, and a class of equivalent scalars with respect to the equation of motion, 
up to the topological scalars of the first class. Therefore in this case also, we can consider a unique 
perfect square that contributes to the dynamics. 

To begin, the more general square made of order 8 scalars has the form : 

{aH{t)‘^ + + 5H{t)H{t) + . 

And its square-root gives the following equation of motion: 

3(a — P + ^ + 25 — 6a) a{tY — 4a(/)d(/)^ 

-I- (7 — (5 -I- 3a) a{tY (3d{t)^ + 4d(/)a^®^(t) -I- 2a{t)a^‘^\t)j = 0. 
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Following the section concerning order 6 , we can say that there are 5 independent contributions 
made of square-roots of order 8 scalars, but in this case, there is also one condition to impose on 
(a, /3, 7 , (5, cr) in order to have second order equations of motion. Therefore there will be only 4 
independent contributions and we can choose them to be : 

^/J 5 = v/-38/Ci - 2448/Cio + 2400/Cii IO 86 /C 12 - 1437Wio - 22QMii + 7927Wi2 SSTWig - 352A4i9 

= 6 v^ {zH{t)H{t) + , 

8 ^ = y/lCi - 144/Cio + 48/Cii -h 27/Ci2 = 24 H{t)^ , 

□i? = = -6 (l2H{tfH{t) + AH{tf + 7H{t)H{t) + ^ 

y/Je = ( 5 /C 1 9(8/Cio - 32/Cii - 7 /C 12 )) = l2VmH{tfH{t), 

where only the last one gives a non vanishing contribution to the equations of motion. We note 
that, of course, because there are much more perfect squares than the previous ones, the fact that 
they form a basis gives once more some non-linear algebraic relations between the order 8 scalars 
for FLRW space-time, and reduces the already reduced FKWC-basis. However it is not the aim of 
this work to find all these relations. We focus here on the unique modified action with respect to 
order 8 scalars : 


^4 




R + - <d(sRf^''R‘^PR^Pf,c.R.pup 







(25) 


Correction to the Friedmann equation 

From this action, we get the acceleration equation : 

+ 2H{t) - 6i^VmH(tf {zH{tf + 4ij(t)) = -87rp. (26) 

In analogy with equations (12) to (14), this last modification leads to an modification of 
the Friedmann equation. Now we are going to see that this unique correction is very interesting 
regarding the problem of the Big Bang singularity as it was shown in Ref m- It is convenient 
to introduce the dimensional quantity pc, critical energy density of the universe, and the constant 
parameter e, such that 

-^ = -2vVm. (27) 

Thus, the modified Friedmann equation induced by the action S '4 is : 

iH{tf + e = 8 ^p. (28) 

The solution that reduces to the standard Friedmann equation in the limit when pc goes to 

infinity is : 

3fl(7=±^('-l + y^ + l), (29) 

Furthermore, if the energy density of the universe is small compared to the critical one, p <C Pc, 
this equation becomes : 

H{tr=^[l-ej-)+Oip^). (30) 

Choosing e = 1, gives the loop quantum cosmology correction to the Friedmann equation m, and 
choosing e = — gives the Randall-Sundrum brane world model [20]. These two corrections are 
regular cosmological solutions and allow to avoid the Big-Bang singularity [Tads]. One should 
also note that the above equation may be obtained within other approches (see mm)- 


(24) 
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4 Static Spherically Symmetric Space-times 


In this section, we study static spherically symmetric space-times, defined by the general metric: 

= -B{r)dt^ + A{r)dr'^ + r'^dn'^. (31) 

First, following exactly the same procedure as for the FLRW case, it is possible to show that there 
is no second order linear combination made of order 6 scalars. The same conclusion is valid for 
the classes 72.2 2 4 - 

Now for order 4 scalars (and more generally for all orders, considering only monomials of the 
curvature tensor), the result of Deser and Ryzhov |18] shows that the most general second order 
action is : 

Sb = j d^x^ [r+ Via ) . (32) 

Moreover, there is no other perfect squares than and B?. Concerning order 4, all 

the scalars that are perfect squares lead to second order equations of motion inside the square. 

Now, starting from order 6 scalars, it is possible to show that there are again only two perfect 
squares. Thus, we can consider the action : 

Sq = j d'^x^ f <5 x/V.RV-i? + ^37 j , (33) 


here is the Cotton tensor, expressed in terms of the Weyl tensor as = —2V‘^iyj^aCT- Its 
square may be written in our basis as = 2 (. 5 f 5 — .ife) — 

In our search for second order differential equations, it is interesting to note that the property 
of order 4 perfect squares is preserved here : these two terms are such that their higher order terms 
cancel perfectly, what makes their associated equations of motion second order. For example the 
term : 


^9 C(r, B{r), A(r),B'{r), A'(r), B"(r), A"(r)) 

— 3r^B{r)^ A{r)A'{r)B"(r) — r^B{r)'^A{r)B'{r)A"{r) + 2r^B{r)^ A{r)^B^^^ (r) |, 


(34) 


where B{r), A{r), B'(r), A'(r), B"{r), A"{r)j is a sum of 15 first order terms (that lead triv¬ 
ially to second order differential equations), is equivalent, up to boundary terms, to the following 
first order expression : 




V3 

6 rB{rYA{rY 


{^B{rfA{rf - AB{rf A{rf - ?>rB{rf A{rf B'{r) 
+ 2r'^B{r)A{r)^B'{r)^ — -r^A{r)‘^B'{r)^ |. 


(35) 


Therefore, we have shown first that, up to order 6 , all the perfect squares that one can build for 
static spherically symmetric space-time are also perfect squares in FLRW as we have seen for the 
case of VctRV'^T?, because in this space-time Wf_i, 9 aa = 0. And secondly, that they also share the 
property that their square-root lead to second order equations of motion for the metric field. 
These perfect squares, coming from the action Sq for A(r) and B(r) are respectively: 

16(7 -f 25)B{rf + 4r( - 67 -f 25)B{rfB'{r) 

+ Ar^{2-i + 5)B{r)B'{rf+r^{--i + 5)B\rf = 0, 

And, 

8(7 -I- 2d) ^ — 1 -I- A{r'^ A{r)B{r)^ + 4 r (57 — 26)B(r)^A'(r) — 18r^^A{r)B{r)B'(r)^ 

+ 8 r (27 -I- S)B{r)^ ^ — rA'(r)B'{r) + A{r) {B'{r) + rB"{r)) ^ (37) 

+ r^{'y- d)B'{r)(5A{r)B'{r)‘^ + 3R(r) (A'(r)R'(r) - 2A{r)B''{r)) ) = 0. 
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The first one provides solutions for B{r), so the two equations decouple. We have found real exact 
vacuum solutions for 3 couples ( 7 , 5) and have computed their associated Kretschmann scalars 
R^vap to see if they suffer from singularites at r = 0 . 

• For 7 = 1 and <5 = 0: 

2 2 
dt^ + dr^ + r^dn^ and ^ 4, (38) 

p + r® 

• For 7 = 0 and <5=1: 

ds^ = -\ df + ^ ^ „ dr'^ + r^dn^ and Rn^cB ^ + ( 39 ) 

1 + pr^ 

• For 7 = 1 and 5 = : 

ds^ = —k dt^ +p dr^ + r^dVt^ and R^''°‘^Rn.ua b — >■ „ V ■ (40) 


where p and k are integration constants. Note that k can be set equal to one by a right choice of 
the time coordinate. 

Recall that the Kretschmann scalar of the Schwarzschild solution of the Einstein equations 
diverges as R^''°‘^ R^ivap —>■ 1/r®. Therefore, we see that our first solution has a milder divergence. 
In the last case, choosing the initial condition p = 1 provides R^''°‘^ R^uap = 0 for all r so there 
is no singularity of curvature because the space-time is flat. However in the second case, as for p 
real, 15p^ — 2p -|- 3 > 0, we cannot ’’cancel” the singularity by a proper choice of initial condition, 
but still, the divergence is again milder than the standard one. 

Now let us note that considering only the two order 8 classes 7^2,2 ’^ 8 . 4 ’ have found a 


unique perfect square 


(that is zero in the FLRW 


2iFiQ — — TVdii “b 2Mi2 — T 4Afi4 — dAIi^ 

case, so possibly again related with Cotton and Weyl ’’conformal” tensors), but its square-root 
does not lead to second order equations of motion, such that we can conjecture that this property 
is true only for the 3 scalars C^'^‘^C^ua and without counting the scalars 

i?^, (D®!?)^, (£ 4 )^, etc, for which the square-root always lead to second order equations. 

To conclude this section, we point out that it would be interesting to find black hole solutions, 
for some values of cr, 7 and 6 , and starting from the following action : 


^7 = J (r + Via ^W^^’'^dW^uaP + \/^ua + ^ ^V.RV^R ) , 


(41) 


where the two last terms would be consider as high energy corrections to the standard Einstein- 
Hilbert action. In this action, the two additional terms are very similar to one present in this S’^. 
We recall that starting from 5 * 5 , one can find exact black hole solutions m- 

After finding black hole solutions for some class of the parameters, it would be possible to 
reduce the number of physically relevant value of (tr, 7 , 6), for example by calculating the Wald 
entropy associated with the action and the class of solutions, like it is done in Ref [23) . by imposing 
the positivity of the entropy to find some contraint on the parameters. 


5 Conclusions 

Considering the order six and eight FKWC-basis of all independent invariant scalars build from 
the metric field and its derivatives, we have found modified gravity models admitting second order 
equations of motion for FLRW and static spherically symmeric spacetimes. 

For FLRW, we have shown that there are unique second order polynomial corrections to the 
Einstein-Hilbert action for both order six and eight. They involve only contractions of curvature 
tensors which is what one could expect to follow from quantum field theory. Yet, concerning static 
spherical symmetry we have seen that there is none. 

Some non-polynomial corrections, that are square-roots of the specific combinations that are 
perfect squares have also been studied for FLRW and spherical symmetry for both order six and 
eight. We have shown that up to order 6 , all the perfect squares lead to second order equations of 
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motion in this way. However, for order eight, there are squares that do not have this property, for 
example the one we have found for spherical symmetry that involves scalars of the TZ 2 2 TZ^ 4 
order eight classes. 

Moreover, concerning FLRW we have seen a ’’mechanism” that provides unique corrections 
with non-vanishing contributions to the equations of motion : there are lots of squares, but their 
square-roots can be decomposed in a very small basis, such that there exist some non-linear alge¬ 
braic relations between order six and order eight scalars. Taken into account, we can choose the 
independent ’’square-root” scalars such that for both order, there are topological scalars and only 
one that is not. 

For order six, the one we choose turned out to be the same as for spherical symmetry. The 
other perfect squares for this last spacetime are squares of the Weyl and Cotton tensors, which 
are identically zeros for FLRW, such that all the squares of static spherical symmetry are present 
in FLRW. Let’s say also that the ressemblance between S'g and S 5 , that admits exact black hole 
solutions, could suggest that it might be possible to find also these kind of solutions from this 
first. It could also be interesting to understand better why the second order corrections specific to 
spherical symmetry involve ’’conformal” tensors. 

We note here that in addition, we have checked if this common non-vanishing square that is 
Vcri?V°’i? is also a square for another physically relevant spacetime, that is Bianchi I, defined by 
the following metric and describing an anisotropic universe as the early universe was [551ISH] : 

-I- a(t)dx^ + I3(t)dy‘^ + ^{t)dz^. (42) 

It turns out that this scalar is also a perfect square in this case and that it leads also to second 
order equations of motions. These results are just copied in the third Appendix, and makes the 
results concerning the action : 

S = j iR + vy^-V^RV^Ry (43) 


our main ones. 

The study of non-polynomial gravity with order eight scalars represents also an interesting 
results in the sense that we can choose the four independent square-root scalars to be three topo¬ 
logical ones, and one that contributes to the equations of motion and gives the well-known 
correction to the Friedmann equation ESIIIT], that leads to a unique physically relevant solution 
for which the limit of small density allows to reproduce the loop quantum cosmology result and the 
one coming from the Randall-Sundrum brane world model. In these two, the big-bang singularity 
is absent and the cosmological solutions are regular ones. 

To finish, we note that our study was done on fixed backgrounds, but General Relativity is 
a background independent theory, so it could be important to check for general background if 
the scalar is still a perfect square, and if it is the case, if its square-root leads in the 

general case to second order equations of motion for the metric field. Moreover, to our knowledge, 
there is no proof that, considering all the FKWC-basis for all orders, it is impossible to find 
linear combinations of invariant scalars that lead to second order partial differential equations 
for a general metric. Indeed the Lovelock papers [131301131] are only restricted their studies to 

lagrangians of the general form: dag^^, dadpg^J^ , i.e. involving curvature tensors, but not 

explicit covariant derivatives of them. Therefore, if this remark is true, one could search for second 
order equations of motion coming from the general lagrangian L(^g^^,dag^i,, ■.■,da...dpg^^. We 
know from our results that there is no such combination for order six and eight scalars because 
there is none for static spherically symmetric spacetime. But as the order grows, even if there are 
more and more terms in the expansions of the scalars, the number of independent scalars inside a 
same class grows as well, such that for some ’’high” order in the FKWC-basis, it might be possible 
to cancel all the higher order derivatives for general metric. 


II 
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A Appendix 

A.l Relations for the reduced FKWC basis 

First we copy the following relations from Ref [23] : 

— 2/li + 8 /I 2 + 24 /I 3 + 2 /I 4 + 24 /I 5 + IBjCs — I 2 /I 7 + /Ig (44) 

is the order 6 Euler density, 

= \r^-2S + T (45) 

is the square of the Weyl tensor, 

Wf = (46) 

and 

Wi = WP-'^^^Wp^p^sW^^P (47) 

are the two independent cubic contractions of the Weyl tensor in four dimensions. 

From the relations of Ref ESI, we have found the following geometrical identities : 

LfivR^'^ = + 4 /I 5 + 4/le — dify + ^■^s) (48) 

Wp^^pS/'^V^RP^ = ifa - ^^2 + + \Cg - (49) 

VPRP'^^^S/pWp^^p = ^7 - 2if5 + 1^8 (50) 

RP^^V^V’^Wp^^p = ^^2 - - ^Le + iis (51) 

= iifs - - ^^8 (52) 

Note that VPV^Lpi, = 0 identically, such that there are no new relations coming from this term. 

Because in four dimensions, £q — 0 and Lp^, = 0, leading to Lp^RP’' — 0 for order 6 scalars, 
the FKWC basis, that does not take into account these relations, is then reduced. 

As for conformally invariant space-times, their metrics verify Wpi,ap = 0 which provides again 
new relations between the scalars coming from : 

RW'^ = 0, Wf = 0, IFf = 0 (note that because £q is a linear combination of Wf and W 2 
[23], the relation £lg = 0 becomes redundant), Wp^ap'^'^^^RP"^ = 0 , VpRP'^°‘^V pWp^ap = 0 , 
RP°'V^V''Wp,jap = 0 and VPWp^ap^'^R''^ = 0. Therefore, for conformally invariant space- 
times, there are 1 relations coming from the corollary of the Lovelock theorem, and 7 coming from 
Wpvap = 0 , so there are 8 scalars less in the reduced basis, i.e. 8 scalars left. 


A.2 General order 6 linear combination for FLRW 

The sum of all order six independent scalars, expressed in terms of a(t) and its derivatives is : 

3 


J = 


{viCi +cci 


a{tf 


(Ji{vi,Xi) a{tY + 'J2{v^,x^) a(t) d{t) d(t) + a^ivi.Xi) a(t) d{t) d{t) 


-I-cr 4 (ui,Xi) a{£f d{t)^ +cr 5 (vi,Xi) a{t)^ d(t)^ +<^6ivi,Xi) a(t)^ d{t) d(t) 

+ a7{vi,Xi) a(t)'‘ a^^^(t)^ -I- as{vi,Xi) a{t)^ d(t)^ + crg{vi,Xi) a{t)'^ d{t) ) 
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with. 


(Ti = 4 :( 2 vi — 2 v3 + 6v4 + 2v5 + 2ve + 6 vy + ISfg + 2x2 + 2:3 + 6x4 
—6x5 — bxQ — 8x7 — 12X8) 

(72 = 2 ( — 2 v 3 + 12 v 4 + 4w5 + 6 x 5 + 24x7 + lOSxs + 30xi + X 2 
— 3 x 3 ~18x4 + 20x5 + 18 x 6 + 32 x 7 + 24 x 8 ) 

(73 = [ — 6 x 2 — 4 x 3 + 24x4 + 14x5 + 6 x 6 + 48x7 + 216x8 + 48xi 
+ 14 x 2 + Txz + 42 x 4 — 20x5 ~ 19 x 6 ~ 86x7 — 12x8] 

< (74 = 8x1 + 2x2 — 8x3 + 24 x 4 + 6x5 + 10x6 + 24 x 7 + 72 x 8 — 12xi — X3 — 6x4 
(75 = —2 ( 18 xi + 5x2 + 2:3 + 6x4 — 4x5 ~ 2x6 ~ 24 x 8 ) 

uq = — ( 36 xi + 8x2 + X3 + 6x4 — 2 x 6 ~ 8x7 + 24 x 8 ) 

(77 = - ( 4 X 5 + 3 X 6 + 4 (X 7 + 3 X 8 )) 

(78 = —2 (6x1 + X2) 

ag = — (12xi + 4X2 + 2^3 + 6x4) 

Now in terms of H{t) and its derivatives, considering only the scalars of the reduced FKWC 
basis : 

J= ^ ViCi+ ^ Xj^j = ‘i\^ai{vi,Xi) H{tf+ a2{vi,Xi) H{tYH{t)+a3{vi,Xi) 
i= 4 , 6,7 4 = 1, 3 , 5 ,8 V 

+ ( 74 (xi,Xi) H{tY + '^b{vi,Xi) +aQ{vi,Xi) H{t)H{t)H{t) 

+ a 7 {vi,Xi)H{if +( 78 (xi,Xi) +a^{vi,Xi) 

with, 

di = 12 (8x4 + 3 (x6 + 4x7)) 

CT2 = 6 ( 24 x 4 + 9 x 6 + 36x7 — 48 xi — 2x3) 

+3 = 4 ( 24 x 4 + 9 x 6 + 30 x 7 ~ 60 xi — 2x3 — 14 x 5 ~ 48 x 8 ) 

+4 = 2 (12x4 + 5 x 6 + 12x7 ~ 24 xi — 2x3) 

< +5 = —7 ( 24 xi + X3) 

CTe = - ( 84 xi + 5x3 + 24 (x5 + 4 x 8 )) 

+7 = -4 (X 5 + 3 X 8 ) 

+8 = - ( 24 xi + X3) 

+9 = - ( 12 xi + X3) 


A.3 VaRV^R in Bianchi I spacetime 

We reproduce only here the value of the particular scalar Vo-T^V'^i? for Bianchi I spacetime, and 
the associated equations of motion of its square-root. 


= -4^a{t)j{t)P(t) |^a(t)/3(t)7(t) + + a(t)i3(t)) j + l3{t) 

+ 'yitf /3{t) - a{t)a{t)i3{t) - a{t){P{t)a{t) + (t))) 

+ /3{tf \a{t)j{t)( a{t)j{t) + a{t)^{t)) + (a{t)a{t) - a{t)a^^'> (t)) + j{t)( aitf'jit) 


— a{t)a{t)^{t) — + a 


n \ 2 


a(f)'‘/3(t)‘‘7(t)^ 


(53) 
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The associated equations of motion, very complicated, yet second order, for respectively a{t), f3{t) 
and j{t) are : 


0 = + a{tyi3(t)j{tyi3{t)l3(t) + a{t)l3{ty-f{ty 




■ 7 (t) {a{typ{t) + a(t)p{t)a(t) + a{t)a{t)j3{t'^ 


+p{ty 


aity^yy+-f{ty(^- 2 a{ty (54) 


■ 3a{t)d{t)a{t)^ + a{tyj{t)j{t)j{t) + a{t)j{ty (^d{tyj{t) + a{t)j{t)d{t) + a{t)d{t)j{t)^ 


0 = -2a{ty^yyp{ty + ayymiyym 




-ayymHty 


+ a{t)j{t)i3y) + d(t) 


+ a{t)l3{t)^{t) 


(55) 


■yyy 


ytfiyf + i{ty[- ^yy + ay)d{t)d{t)j + ayy-iy)^y)^y) 


0 = -ayy-fitypyy + ayy i3y)-iyy pyyyy.) + pyy 


Mtfiyy 


+ lyy (“ + a{t)dy)dy,)^ 

+ otyyiy)iy){ay)^y) + 3a{t)^{t)j + ayy-iyy(^^y)ay) + dy)^y)j 


■ ayypyy^it) 


aiyyyyyf+iy) 


■ ay)jy)i3y) + a 


(56) 
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